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Using the DMRG method we calculate the surface spin correlation function, Cl{1) = {SI Sf^^i_i) , 
in the spin 5* = 3/2 antiferromagnetic Heisenberg chain. For comparison we also investigate the 
5* = 1/2 chain with 5 = 1 impurity end spins and the 5* = 1 chain. In the half-integer spin models 
the end-to-end correlations are found to decay to zero logarithmically, Cl(1) ~ (logL)"^"^, with 
d — 0.13(2). We find no surface order, in clear contrast with the behavior of the S — 1 chain, where 
exponentially localized end spins induce finite surface correlations. The lack of surface order implies 
that end spins do not exist in the strict sense. However, the system possesses a logarithmically 
weakly delocalizing boundary excitation, which, for any chain lengths attainable numerically or 
even experimentally, creates the illusion of an end spin. This mode is responsible for the first gap, 
which vanishes asymptotically as Ai « {-Kvsd) /{L\n L), where vs is the sound velocity and d is the 
logarithmic decay exponent. For the half-integer spin models our results on the surface correlations 
and on the first gap support universality. Those for the second gap are less conclusive, due to strong 
higher-order corrections. 

PACS numbers: 75.10.Jm, 75.40.Cx, 68.35.Rh 
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I. INTRODUCTION 

One of the most studied models of interacting quantum 
systems is the antiferromagnetic Heisenberg chain. It is 
defined by the Hamiltonian 



L-l 



(1) 



where Si is the spin-S* operator at site i. It is known since 
Haldane's seminal worki that the low-energy properties 
of the model for half-integer and integer values of the 
spin are different. Half-integer spin chains are gapless, 
and the bulk correlations are quasi-long-ranged. On the 
other hand, integer spin chains have a finite gap, there 
is a hidden topological order yS. and the bulk correlations 
are short-ranged, involving a finite correlation length ^. 

The universality class of half-integer spin chains was 
debated for some time, by now it is generally accepted 
that all are in the same bulk universality class, indepen- 
dent of the value of the spini^i^ A thorough numerical 
investigation for the 5 = 3/2 chain& verified convincingly 
that the bulk spin-spin correlation function {SqS^) ^ r~'' 
decays with the same bulk exponent 77 = 1 as that of the 
5 = 1/2 model. It is a priori not clear, however, to what 
extent this universality can hold for the surface behavior. 

In an open chain with L sites it is of interest to measure 
the two-point function, 



Cl{1) ~ (5f52+i_() 



(2) 



where the two spins are positioned symmetrically to the 
middle of the chain. When 1 <C ^ < i/2 the two spins 
are far from the boundary and the scaling behavior is 



determined by the bulk exponents. For instance, Cl{olL) 
with < a < 1/2 scales as CL{aL) ~ i"''. 

On the other hand, for I ^ L fixed, the two spins 
are close to the boundary, and Cl{1) probes the surface 
critical behavior. In particular, for / = 1 we obtain the 
end-to-end correlation function Cl(1). For the 5 = 1/2 
chain, as obtained by analytical methods, the end-to-end 
correlation is also quasi-long-ranged, but the asymptotic 
decay, Cl(1) ~ L~^% involves the so-called surface ex- 
ponent rjg = which differs from the bulk exponent 

For higher 5 there may be an even more drastic differ- 
ence between bulk and surface correlations; the surface 
correlations may become long-ranged. This phenomenon 
is best understood for the 5 = 1 chain, in which there is a 
well-established surface order, C(l) = limi^oo C'l(I) > 
0. In general the order decays exponentially with the 
distance from the surface. 



C{1) = lim Cl(0 -exp(-2Ve.) • 



(3) 



where is the surface correlation length, which is equal 
to the bulk correlation length ^ of the 5=1 chain&S 
(Note, however, the factor 2 in the exponent, which is 
present due to the fact that a unit change in I takes the 
two spins two sites closer to each other.) The leading 
finite-size correction to the correlation function is also 
exponential; 



CL{l)~C{l)^eM-L/if) ■ 



(4) 



Here again the length characterizing the finite-size cor- 
rections is expected to be equal to ^, thus there is only 
one length scale in the system. Note that we have ne- 
glected possible algebraic prefactors in Eqs. (j3Jl and Q). 
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The origin of surface long-range order in the S = 1 
chain can be understood in the valence-bond-solid (VBS) 
framework^ Representing each S = 1 spin by two 
symmetrically coupled spin-1/2 degrees of freedom, the 
ground state can be written as the product of singlet 
bonds connecting neighboring spins, except for two re- 
maining localized spin-l/2s at the ends. These localized 
end spins are responsible for the behavior of surface spin 
correlations. The VBS construction is exact for the S ~ 1 
Affleck-Kenedi-Lieb-Tasaki (AKLT) model. 



L-l 

Haklt = ^ Si ■ Si+i 
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(5) 



However, numerical studies have confirmed the above sce- 
nario for the pure 5 = 1 chain. Hi, tooiS 

In order to understand the behavior of open spin chains 
with 5* > 1, Tai-Kai Ng proposedii a theory based on 
an extension of the VBS picture. According to this, each 
spin-S' variable can be replaced by 25 symmetrically cou- 
pled spin-l/2s, among which singlet bonds are formed. 
For integer S all bulk spins are involved in singlet for- 
mation, except for those at the chain ends, which form 
localized edge spins with spin 5/2. Surface spin correla- 
tions are then expected to have the same characteristics 
as for 5 = 1, summarized in Eqs. l|3Jl and @. The ground 
state becomes (5 -I- 1)^ times degenerate in the thermo- 
dynamic limit. For finite chains this degeneracy is lifted 
by an exponentially small gap, which is due to the ex- 
ponentially small effective interaction between the end 
spins. Above the degenerate ground states the spectrum 
is gapped (Haldane gap), and the finite-size corrections to 
these levels are algebraic (note the boundary condition) 
DMRG-studies of the 5 = 2 spin chain have shown evi- 
dence in favor of the above scenario, although the corre- 
lation length is much longer than for 5 = 

For half-integer spins, this scenario implies that there is 
one spin-1/2 at each bulk site which is not involved in sin- 
glet formation. Moreover, there remain spin-(5 -I- l/2)/2 
degrees of freedom at the chain ends. Consequently, as 
Ng arguesjii the effective low-energy model is a spin-1/2 
chain (a Luttinger hquid) perturbed by (5 -I- l/2)/2 im- 
purity spins coupled antiferromagnetically at the ends: 



(S+l/2)/2 
/2 



(A>0), 5 =half-integer, (6) 



where 



L-2 



K/2(A) ^Y.^^- ^'+1 + A • 52 + 5l_i • ai) , (7) 



with 5 ((j) being spin-1/2 (spin-cr) operators, and A > 
is some unknown, but positive, impurity coupling. 

The 7i^y2 quantum impurity problem was analyzed 
earlier by Eggert and Affleckii^ They argued that the 
surface spin is partially screened by the Luttinger liquid, 
leading to an effective end spin of value aend = cr — 1/2. 



This remaining spin is coupled ferromagnetically to the 
bulk, thus 



K/2(A>0) 



nlj,"\~g < 0), 



(8) 



as far as the low-energy behavior is considered. How- 
ever, the coupling g — g{L) > is marginally irrelevant, 
satisfying a Kondo-type RG equation 



dg 
d\nL 



-9 



whose solution is 



9iL) 
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In(BL) 



(9) 



(10) 



Here B is an appropriate inverse length depending on the 
bare coupling. Note that B can be rather different from 
unity, thus its effect cannot be neglected for moderately 
long chains. 

Equation H1U|) implies that end spins become asymp- 
totically free, g{L) ^ for i ^ oo. In the thermody- 
namic limit the spectrum is that of an 5 = 1/2 chain with 
L — 2 sites and two decoupled spin-CTend spins. This gives 
rise to an (2o'e„£;+l)^-fold quasi-degeneracy of the ground 
state. For finite length the quasi-degenerate states are 
separated from the ground state by an energy 



Adcg(i) 
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L ^ L\n[BL)' 



(11) 



The other higher-lying excitations in the spectrum are 
predicted to be delocalized excitations of the Luttinger 
liquid. 

Some of these field-theoretical predictions have been 
checked numerically for half-integer spin Hcisenberg 
chains ii2iiiiSiiLi& In particular, indications for the pre- 
dicted number of edge states were found, and the corre- 
sponding excitation energies were shown to vanish faster 
than 1/L. Higher- lying states were demonstrated to scale 
as 1/L. However, several aspects of the system have re- 
mained unclear. The accurate form of the excitation en- 
ergy of the edge state lacks an accurate analysis. It is 
not clear whether the logarithmic term predicted by the 
quantum impurity calculation is indeed consistent with 
the numerical results. Also, in order to test universal- 
ity the structure of the higher lying excitations requires 
additional analysis. 

Another open question involves the surface correlation 
function, and in particular the surface long-range order. 
Indeed, the VBS analysis is based on the assumption that 
the bulk chain possesses a short-range RVB "backbone" 
which remains inert for the low-energy behavior. Its only 
effect is to provide localized edge spins. This picture 
would imply that this backbone mediates correlation be- 
tween the two ends of the chain in much the same way 
as it does for the 5 = 1 Heisenberg chain. However, 
the presence of the extra spin-1/2 component in the bulk 
may have an important effect, and it is a priori unclear 
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whether this effect is weak enough not to destroy the 
surface order. In Ng's theory the edge spins are local- 
ized to the ends, and their effective coupling scales to 
zero faster than 1/L. When the chain length is even, 
the lowest-lying state is a singlet, in which both the bulk 
of the chain and the two end spins form singlets inde- 
pendently. Thus there is a clear anticorrclation between 
end spins. Localization of the end spins means that they 
contribute a finite amount to the first and last spins of 
the chain, thus this anticorrclation assures finite surface 
order, C{1) > 0. 

Our aim in this paper is to elaborate on these questions 
in the case of the S' = 3/2 antiferromagnetic Heisenberg 
chain. We study the spin-3/2 chain, 7^3/2, using the 
DMRG method, and calculate low-lying excitation en- 
ergies and surface correlation function. For comparison 
with the VBS theory, we also consider the S = \ chains, 
Ti-AKLT and Hi, and the effective model, Hy2i^ > 
which is predicted by Ng's theory to belong to the same 
universality class as the 5 = 3/2 chain. 

The structure of the paper is the following. In Sec. UTI 
we discuss the high precision DMRG method used in the 
calculations. The numerically calculated excitation en- 
ergies and the surface correlation function are analyzed 
in Sec. IIIII and in Sec. IIVI respectively. A discussion of 
the results is presented in Sec. There is an Appendix 
attached to the paper, where we list the exact surface 
critical properties of an Ising quantum chain with an ex- 
tended defect, which serves as an analogy with rigorous 
results. 



II. NUMERICAL METHOD 

In the calculations we used open chains of various fi- 
nite lengths, L = 12, 24, 36, ...,216, and calculated en- 
ergy gaps and the surface correlation function using the 
DMRG method (for recent reviews sec Ref. ITfll) . We per- 
formed the DMRG calculations by using both the stan- 
dard technique, i.e., keeping the number of block states 
fixed, and the dynamic block state selection (DBSS) 
approach All eigenstates of the model were targeted 
independently using two or three DMRG sweeps. 

In the standard procedure 500 — 800 block states were 
used. It was found that for the largest systems {L = 192 
for the S = 3/2 chain and L = 216 for the S = 1/2 
chain with impurities), the truncation error varied in the 
range lO"'^ - 10"**. The DBSS approachSa allows for 
a more rigorous control of numerical accuracy, and we 
set the threshold value of the quantum information loss 
X to 10~^. The minimum number of block states was 
set to 256. The entropy sum rule was checked for all 
finite chain lengths for each DMRG sweeps, and it was 
found that the sum rule was satisfied after the second 
sweep already. The maximum number of block states 
varied in the range 600 — 1400. In order to have high 
precision results, we set the residual error of the Davidson 
diagonalization algorithm to 10~^. 



III. EXCITATION ENERGIES 

The spectrum of the 5 = 1 chain with open boundaries 
is well known in the literaturei^^ The first gap (singlet- 
triplet gap) , which is related to the localized surface exci- 
tation, is exponentially small in L. In contrast with this, 
the second gap (singlet-quintuplet gap), which describes 
a bulk excitation, is finite and this defines the well-known 
Haldane gap. 

For S = 3/2, numerical results on the spectrum of open 
chains were presented in Ref. ll2L in which the authors 
aimed to make a comparison with the field-theoretical 
prediction of Ngjii Indeed the numerical results of Ref. 
I12I are compatible with the prediction that the first gap 
vanishes faster than 1/L. However, as their quantitative 
analysis of the finite-size dependence of the first gap was 
rather sketchy, we are going to repeat this type of anal- 
ysis here. Furthermore, we make a comparative study 
of the spectrum of the effective model 7il^2W' which 
arises from the VBS picture. For A, which is an a priori 
unknown positive bare coupling, we test three different 
values, A = 0.5, 1.0 and 4.0. We expect that the actual 
value does not matter due to universality. 

In both half-integer spin models the remaining effec- 
tive end spins are aend = 1/2, thus the quasi-degeneracy 
involves the singlet ground state and a triplet state. Thus 
the first gap, defined by the lowest state in the S^^f = 1 
sector, is expected to take the form of Eq. Hll|) . The 
higher-lying excitations are delocalized excitations. They 
scale as 1/L and form conformal towers. In the follow- 
ing we analyze the gaps A„, n = 2, 3, . . . , defined by the 
lowest states in the sectors 5^^^ = n. 

Conformal invariance dictates that the gap A„, asso- 
ciated with the scaling field (j)n, satisfies 



A„ = 



TTVS 



dn 



In(BL) 



O 



InlnL 
In^L 



(12) 



where vs is the non-universal sound velocity and Xn is the 
scaling dimension of (/)„, defined as {(f>n{0)4'n{r) ~ r~^^". 
The sound velocity is known exactly for the 5 = 1/2 
chain: Vi/2 = Tr/2m^ For 5 = 3/2 the best numerical 
estimate is W3/2 = 3.87& 

The logarithmic correction appears due to the pres- 
ence of marginal operators 0j„ai.gi24 The coefficient dn 
is universal, determined unambiguously by the operator 
product expansion (OPE) rules of the underlying field 
theoryi^^iS^ Assuming that the OPE reads schematically 
as 



0n • (t^n = bn 0marg + • ■ • , (13) 

the coefficient of the log correction turns out to be24 

2bn 
dn = —[-■ 

b 

Other, more irrelevant operators may also give correc- 
tions but these decay faster than the leading order loga- 
rithmic correction. In the present case there are, in fact. 



(14) 
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two marginal operators in the theory: the boundary per- 
turbation represented by the effective edge spin discussed 
above, and the standard bulk marginal operator present 
in all half-integer spin chains In the following we assume 
that even in this case Eq. (|12|l remains valid, with dn rep- 
resenting an effective value, arising due to the combined 
effect of the two marginal operators. 

Note that the expected form of the gaps in Eqs. 1)11(1 
and (|12|) are formally the same. We can think about the 
first gap as satisfying Eq. H12() with scaling dimension 
xi — and an appropriate constant di . 

Neither for 7^3/2 nor for the Tij^yjl-^) model are we 
aware of any rigorous analytical predications for Xn and 
dn- Nonetheless, if the two models are in the same (sur- 
face) universality class, the full operator content should 
be identical. A possibility, suggested implicitly in Ref.H^ 
is that the critical exponents associated with the higher- 
lying excitations in the spectrum arc the same as the ex- 
ponents in the pure 3=1/2 model with open boundary 
condition, Ti.i/2- If this were true, the edge spins would 
only effect the lowest gap Ai but not any higher-lying 
excitations. For the 7ii/2 chain the operator content is 
rigorously known 



dn = -n{n + l)/2. 



(15) 



However, numerical results of Ref.ll2l when adjusted cor- 
rectly with the sound velocity U3/2J do not support con- 
vincingly this scenario for the S" = 3/2 chain. It may well 
be that the end perturbation, which has similar effect to 
a modified boundary condition, substantially alters the 
surface operator content. This phenomenon is in fact 
possible: Ti.i/2 (open) and TI1/2 + Sl ■ Si (periodic) have 
different operator contenti£iS& Nevertheless, we can use 
the values in Eq. H15|) as reference values for comparison 
in the sequel. 



A. Analysis of the first gap 

The size dependence of the first gap is shown in Fig. 
n both for 7^3/2 and for the effective model Hy2W 
with three different values of the impurity coupling, 
A = 0.5, 1.0,4.0. The data seemingly fit onto a straight 
line in the log- log plot, with an average slope larger than 
1. To obtain a more quantitative description, we have 
calculated the local slopes from two point fits. This cor- 
responds to the effective gap exponent C{L), defined as 
Ai(L) ~ L~''(^\ and the true gap exponent is obtained 
as C = limi^oo C(-^)- 

For the S = 3/2 chain, the effective exponent first in- 
creases with L, then passing through a maximum around 
L = 60, it starts to decrease rapidly. The same data for 
the impurity models are monotonically decreasing with 
L. To analyze the size dependence of the effective gap 
exponents we use the field-theoretical result in Eq. (|ll|l . 
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FIG. 1: Log-log plot of the first gap Ai vs. L; The inset 
shows the associated local exponent C(^) vs. l/ln(i3L) with 
B chosen appropriately for each data set. 
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FIG. 2: Length-dependent effective parameters, -B(I/) and 
di{L) determined from two-point fits. Dotted line segments 
indicate uncertainty. 



which predicts logarithmic corrections. Therefore the ef- 
fective exponents are presented in the inset of Fig. da) 
as a function of l/ln(i?ij) with the optimal B chosen in- 
dependently for each data set (see later). For almost all 
curves the behavior seems to be nicely compatible with 
the expected result C = 1 with logarithmic corrections, 
except for A ~ 0.5 where the asymptotic regime has prob- 
ably not yet been reached. 

Assuming that the first gap indeed scales according 
to Eq. I|12|) with xi = 0, the two remaining parame- 
ters B and di can be determined from the numerical 
data. The consistency of these parameters with other 
results can serve as a justification for the fitting ansatz. 
In Fig. [3 we present the effective, size-dependent param- 
eters B{L) and di{L), determined from two-point fits. 
The parameter B which depends on the bare couplings 
of the marginal operators is not universal. Indeed we see 
that it varies considerably for the four models we con- 
sider. Our best estimates for B = limi^oo B{L) are as 
follows: 



0.11 ±0.02 iorH 



B 



0.25 ±0.1 
1.2 ±1 

4±2 
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Note, however, that B is an extremely sensitive param- 
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FIG. 3: Scaled first gap LAi as a function of l/ln{BL), witii 
B cliosen individually for each curve. The common slope 
(dotted line) determines the universal value di = 0.13 ± 0.02. 



FIG. 4: Scaling dimension X2 associated with the scaled sec- 
ond gap A2L as a function of 1/ In BL. The dotted line indi- 
cates the conjectured scaling. 



eter since it is in the argument of a logarithm. Even 
the 10"^" relative precision of the gap values may not be 
enough for a reliable extrapolation. Indeed, we cannot 
exclude that the sharp downturn of B(L) for the impu- 
rity models above L « 72 may indicate a systematic error 
stemming from the loss of numerical precision, and not 
a real physical effect. Thus the determination of B for 
these models bears a rather large error. 

Even though the error of B is considerable, many other 
quantities do not inherit this large error bar due to their 
logarithmic dependence on B. This is clearly seen for di 
in Fig.El Our best estimate for H3/2 is di = 0.14 ±0.02, 
and for TiJ^j "^i = ^.12 ± 0.02. The deviation is small, 
and a common value at di = 0.13 ± 0.02 is consistent 
with our error estimates. This supports the conjecture 
that these models belong to the same universality class. 
Note that di is positive, whereas for the pure 5* = 1/2 
chain, 'H1/2, d„ in Eq. H15|l is negative for all n. 

Using the values of B determined in Eq. (|16|) . the 
first scaled gap is convincingly linear as a function of 
l/hi{BL) as shown in Fig.|2| The common slope is given 
by di within the error bar. 

The corresponding value di for the 5 = 3/2 chain can 
also be calculated from the results of Lou et ali^ on 
the coefficient of the log correction. Their data yields 
di = 0.32, however, they used a very rough fitting with 
B = 1, which is very far from our estimated value. As for 
B itself, our value can be compared to the one obtained 
in Rcf. iQ, namely B = 0.6 ± 0.1. The two values differ 
by a factor of five. Notice, however, that Ref. |^ used 
periodic boundary condition, and thus measured the log- 
arithmic corrections due to the bulk marginal operator 
only, whereas our B is expected to represent the com- 
bined effect of both marginal operators. Also B{L) is 
monotonically decreasing for 5* = 3/2, thus determining 
it from shorter chains such as L = 60 used in Ref. |^ may 
yield a higher effective value. 



B. Analysis of the second gap 

Our numerical data indicate both for the S = 3/2 
chain and for the impurity models that the second gap 
A2 scales as 1/i, and the associated dimension X2 is fi- 
nite. The scaled second gap L/S.2/ {ttvs) is shown in Fig. 
2] as a function of l/\ii[BL). For B we have used the 
values determined from the first gap. It is clear that the 
finite size corrections are very strong and the asymptotic 
regimes have not yet been reached. 

The curves for T-l\^^{X) go through a local minimum, 
then increase rapidly as L increases. A limiting value 
around 12 ~ 1 seems plausible, but the approximation 
error is large, and in fact nothing above 0.8 can be firmly 
excluded. As mentioned above the scaling dimension 
X2 = 1 would be expected if the delocalized excitations 
behaved as those of the pure H1/2 model for which the 
lowest scaling dimension given by Eq. (|15|l is 1. 

The coefficient of the leading logarithmic correction 
d2 may again be universal in accordance with Eq. H12|l . 
The most probable value, consistent with our data is 
d2 ~ ~2.5±1. This has different sign than di. It also dif- 
fers considerably from di = — 1 characterizing the lowest 
gap in 7^1/2. Nevertheless, we should warn the Reader 
that the estimation of X2 and d2 depends very much on 
the precision of B. Changing B practically acts as a hor- 
izontal shift in the Figure which can modify substantially 
the extrapolation results. 

For the 5 = 3/2 chain, as seen in Fig. 21 the finite- 
size dependence is very small. The scaled gap initially 
decreases in L, then starts to increase very slightly for 
larger L. If this were the thermodynamic limit, we would 
conclude that X2 ~ 0.47 ± 0.01. This would contradict 
universality with the impurity model. However, having 
seen the nonmonotonic behavior for the impurity model, 
there is some doubt for this conclusion. We have to leave 
open the possibility that the local minimum of X2 may 
be followed by a strong upturn for larger values of L. 
We estimate that the crossover length for this upturn 
should be, in worst case, around l/ln(BL) w 0.2, i.e., 
L ~ 1000, which is just slightly beyond the system sizes 
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we can currently attain with reasonable accuracy. 



IV. SURFACE CORRELATION FUNCTIONS 

In this Section we investigate the surface correlation 
function Cl{1)- Wc start our analysis with the S — 1 
chains, TYaklt and Tii, where the surface order and the 
localization of edge spins are well-understood. These 
models serve as benchmarks, and we ask to what extent 
these results can be generalized to the half-integer spin 
models, ^3/2 and 7iJ^2('^)- 



1. S = 1 spin chain 

Correlation functions can be calculated rigorously for 
the AKLT model, Haklt, where the exact ground state 
is the VBS state iSiSLSSiSa It turns out that the two point 
correlation function between any two points i and j > i 
in a finite open chain with length L in the singlet ground 
state is 



(S'f 5*1) AKLT 



(17) 



3[l-(-i)--i] 



L-j+i 



From this the surface correlation function with i ~ I, 
j = L + 1 — I reads 



Cl(0 



(18) 



3[i-(-i)M 



L + l-2i 



2/-1' 



For L —> 00 there is a finite surface order 
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C{1) = -4 Q j = -4e~2i/«AKL 



(19) 



with ^AKLT = 1/ ln(3) the (bulk) correlation length of the 
AKLT model. Note that C(l) = Wuil^^ Cl(1) = -4/9, 
which is exactly equal in magnitude to the string order 
parameter. This is a peculiarity of the AKLT model, 
but there is a general qualitative relationship between 
the existence of end spins, the string order in the bulk, 
the surface order, and the non-vanishing entanglement 
between the end spins all holding for the pure S = 
1 chain too. All these quantities can be used as order 
parameters in the Haldane phase. 

As expected, the /-dependence of C; and also the finite- 
size corrections are determined by the bulk correlation 
length ^AKLT- For instance, for Ci(l) we find the expan- 
sion 



Cl{1) 



L^p — Lj Jaklt 
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(20) 
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FIG. 5: Exponential finite-size scaling of the surface correla- 
tion function Cl{1) of the S = 1 HAF chain for I = 1,2,3. 
Inset shows the behavior of the extrapolated values C{V) as a 
function of I. 



For Til we determined the surface order numerically. 
Precise numerical estimates for C(Z) and ^/ are obtained 
in the following way. Assuming finite-size corrections in 
the form of Eq. we obtained effective, finite-size esti- 
mates from three-point fit both for C{l,L) and £,f{L). 

The fitted values for C(l, L) converge very rapidly, and 
we find a clear surface order C{1) = liniL^oo C{1, L) > 0, 
whose value for / = 1 reads C(l) = -0.28306484(1). 
Note for comparison that the string order parameter 
in this case is 0.374335iS The correlation length has 
a stronger L dependence, which can be a result of a 
neglected algebraic prefactor in the r.h.s. of Eq. |0J. 
Its extrapolated value for the end-to-end correlation is 
^/ = 6.1(1), which should be equal to the bulk correla- 
tion length ^ = 6.03|^ as expected on general grounds. 

The finite-size corrections to the surface spin corre- 
lations are illustrated in Fig. O in which Cl{1) — C{1) 
for I = 1,2,3 is plotted function of L in a semi- 
logarithmic scale. Having an appropriate extrapolation 
value C{1), the points lie on an approximately straight 
line, the asymptotic slope of which defines the inverse of 
The slope of all lines in Fig. ISjare compatible with 
the bulk correlation length. 

The extrapolated values of the surface correlation func- 
tion C(Z), obtained for different l-s are plotted in the in- 
set of Fig. Isjin a log-log plot. The asymptotic slope is 
again compatible with the known bulk correlation length. 
An accurate estimate is obtained by calculating effective 
^(Z) values by two-point fits, see Table HJ which is then 
extrapolated for / cxo as = 6.05(3). 



2. S* = 3/2 spin chain 

The numerically determined values of the surface cor- 
relation function Cl{1) for I = 1,2,3 are given in Table 
nil We believe that these values are accurate up to at 
least 8 digits. 

The fundamental question we would like to answer is 
whether C{1) — linii^oo C'l(/) remains finite as for the 
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TABLE I: The surface correlation function, C{1), of the S — 
1 chain and the (surface) correlation length, obtained by 
two-point fit. The extrapolated values, are obtained 

by the BST convergence acceleration algorithm—*— with the 
parameter, to — 2.5. 





C{1) 




Ss 


1 


0.28306484 






3 


0.13931623 


5.64230917 




5 


0.06887108 


5.67770491 


5.69806489 


7 


0.03439726 


5.76153408 


5.84348388 


9 


0.01734266 


5.84105901 


5.95103010 


11 


0.00880481 


5.90082421 


6.00726945 


13 


0.00449158 


5.94270505 


6.03378854 



CJ3) ' 




1 

■ / / 




11 / 




I/CJ2) 





0.02 0.04 0.06 




0.02 0.04 0.06 



FIG. 6: Local fit parameters for Eq. 12111 calculated from 
three-point fits. 



TABLE IL The surface correlation function C'l{1) of the S — 
3/2 chain for / = 1,2,3. 



L 


Cl{1) 


Cl(2) 


Cl(3) 


12 


-0.6149116033 


-0.3463729630 


-0.4577771615 


24 


-0.5312680335 


-0.2796013080 


-0.3483003563 


36 


-0.4960012887 


-0.2563400653 


-0.3154186768 


48 


-0.4754275850 


-0.2437536782 


-0.2984765312 


60 


-0.4614763573 


-0.2355733053 


-0.2877411795 


72 


-0.4511552246 


-0.2296862928 


-0.2801355963 


84 


-0.4430747706 


-0.2251658062 


-0.2743575096 


96 


-0.4364949818 


-0.2215368743 


-0.2697544487 


108 


-0.4309814145 


-0.2185291971 


-0.2659614093 


120 


-0.4262589357 


-0.2159752350 


-0.2627549925 


132 


-0.4221438229 


-0.2137653001 


-0.2599904263 


144 


-0.4185068985 


-0.2118233164 


-0.2575681425 


156 


-0.4152559777 


-0.2100956569 


-0.2554183252 


168 


-0.4123218910 


-0.2085426652 


-0.2534897645 


180 


-0.4096518527 


-0.2071345489 


-0.2517442428 


192 


-0.4072038977 


-0.2058474800 


-0.2501511731 



S = 1 chain or scales to zero. Since bulk correlations 
are algebraic, we expect that neither Cl{1) can contain 
a finite characteristic length, in clear contrast to the be- 
havior of the 5=1 chain. 

To obtain possible working hypotheses on Cl{1), let us 
briefly consider an exactly soluble model, the Ising quan- 
tum chain with an extended defect^ whose properties are 
summarized in the Appendix. For this model there are 
three different type of surface critical behavior: 

i) For weak enough defects C(l) is zero and the end- 
to-end correlations vanish algebraically. In this case the 
scaled gap, LAi, has a finite limiting value predicted by 
conformal invariance. 

ii) For a critical strength of the defect C{1) is zero and 
the end-to-end correlations decay logarithmically. In this 
case the scaled gap vanishes logarithmically. 

iii) Finally, for strong enough defects C{1) is finite, 
which is approached by algebraic corrections. At the 
same time the scaled first gap goes to zero algebraically 



with L. 

We start to analyse surface correlations for the S = 3/2 
spin chain using an algebraic ansatz 

Cl(0 , (21) 

which works for regions i) and iii) of the inhomogeneous 
Ising quantum chain. Local approximations to the three 
fit parameters can be determined from three-point fits us- 
ing the consecutive numerical values for Cl-i2,Cl and 
Cl+12- These L-dependent local parameters are depicted 
in Fig. El as a function of 1/L. If the ansatz were good, 
the fit parameters should stabilize as L increases. How- 
ever, as is clearly seen in the Figures, all parameters are 
monotonically decreasing functions of L, and there is no 
sign of convergence. The effective exponent oi{L) reaches 
0.3 for L = 192, and it continues to fall rapidly. The 
effective values of a are much smaller than the known 
surface exponent of the S ~ 1/2 chain. In fact, the 
rapid decrease seems to be consistent with a vanishing 
exponent in the L — > oo limit. Approximants for C(l) 
also decrease considerably for increasing L. The most 
probable interpretation is that the ansatz in Eq. (|21|l is 
inadequate. Our experience was similar with Cl(2) and 
Cl(3) as is shown in Fig. El 

Now we try to analyze the data instead of Eq. H21|) 
with a purely logarithmic behavior with vanishing C{1), 

CUn^-AilniBLy^, (22) 

which contains two free parameters Ai and fi, if we use 
the value of B determined from the first gap. This form 
corresponds to the borderline case ii) for the inhomo- 
geneous Ising quantum chain. This type of analysis is 
supported by the fact that the scaled first gap is van- 
ishing logarithmically in both models. Indeed, this fit 
describes the correlation data more satisfactorily, as is 
obvious from the stability of the local parameters in Fig. 
□ 

In order to present a theoretical basis to this procedure 
we note that having a marginal surface operator which 
scales as in Eq. H10|l . the end-to-end Green's function, 
^l"'' ~ {4>n{^)4>n{L)) ^ associated with operator (/>„, is 



0.02 0.04 0.06 0.08 0.1 
1/L 



0.02 0.04 0.06 0.08 0.1 
1/L 



0.02 0.04 0.06 



0.02 0.04 0.06 



FIG. 7: Local fit parameters for Eq. 11221 calculated from two- 
point fits for the 7^3/2 model. 



believed to satisfy the Callan-Synmnzik equatio; 



5.34 



Gi") = 0, 



(23) 



where the beta and gamma functions are 

f3{g) =: -TTbg^; 7„(5) = x„ + 27r6„5r. (24) 
From this, up to 0{g), the Greens function reads^ 



G^"^ = const L-^-" 



1 



-1 -46„/6 



nbg{L) 



(25) 



We have seen that for the first excited (boundary) state 
xi = 0. Using then the actual scahng form of g{L) in Eq. 
(|10|l . we obtain 



G«^ln(i?L) 



(26) 



where we have used that by Eq. ^fT^ di = 2bi/bA Ex- 
pecting that for smaU I the leading contribution to Cl{1) 
is given by G^^'' we obtain the fitting ansatz in Eq. (|22|l . 

The results for Ai and ^ are shown in Fig. [7| For 
I = 1,2,3 the measured exponent fj, is estimated to be 
ji = 0.31 ± 0.02. This value is convincingly close to 
2di = 0.28 ± 0.04, obtained earlier from the gap, sup- 
porting the theory outlined above. We note that the 
exact results of the inhomogeneous Ising quantum chain 
in case ii) are in complete agreement with the predicted 
behavior. The prefactor of the logarithmically vanishing 
scaled gap in Eq. ljA5|l is just the half of the exponent 
of the logarithmic decay of the end-to-end correlations in 
Eq. ^K^. 



3. S = 1/2 spin chain with impurity end spins 

Finally we analyze the surface correlations for the im- 
purity model T-L\i2{^) with A = 1. The fit parameters for 



FIG. 8: Local fit parameters for Eq. 1221 calculated from two- 
point fits for the T-L\j2{^ ~ 1) model. 



Eq. H22() are presented in Fig.|Sl Scaling to finite values is 
very convincingly satisfied, especially for the end-to-end 
correlation, Gl(1). There the extrapolated exponent is 
jjL = 0.21 and the amplitude is A = 0.54. The exponent is 
found to be the same with somewhat higher error bars for 
any I in Cl{1), and the amplitude Ai is a monotonically 
decreasing function of /. Again the exponent nicely sat- 
isfies the expected relation /i = 2di with di = 0.12 ±0.02 
from the first gap. 



V. DISCUSSION 

In this paper we have studied the energy spectrum and 
the surface spin correlations of the S* = 3/2 antiferromag- 
netic Heisenberg chain. The analysis of the numerical 
data for both quantities was difficult due to strong loga- 
rithmic terms, which are the consequence of the presence 
of a bulk and a surface marginal operator in the fixed 
point. Our very accurate numerical data calculated by an 
advanced DMRG algorithm on finite open chains made 
it possible to infer the true asymptotic behavior. The 
end-to-end correlations are found to approach a vanish- 
ing limiting value, however the decay is extremely slow, 
it is logarithmic in L. This type of logarithmic decay is 
accompanied by a logarithmically vanishing scaled first 
gap in the spectrum of the system. This latter result is in 
accordance with a field-theoretical analysis of the prob- 
lem, which is based on the assumption of localized edge 
states as in the 5* = 1 chain. Localization of edge spins, 
however, would imply an exponentially decaying surface 
correlation function, which is in contradiction with our 
numerical results. We can thus say that the edge states 
are logarithmically delocalizing, or framing it differently, 
are just "quasi-localized" . Having in mind the analogous 
result for the inhomogeneous Ising quantum chain in Eq. 
(|A6p we can say that the edge spins are concentrated 
in a region of size InL at the end of the chain. The 
delocalization of edge spins is so weak that even for a 
macroscopic system we may have the illusion that they 
exist and can be measured experimentally. For instance, 
we have seen that for L = 108 the end-to-end corre- 
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lation is C'l(I) « —0.43, while for a macroscopic sys- 
tem with L = 10^^ spins the fitted scaHng form imphes 
Cl(1) w —0.18, still a relatively large value. 

For comparison we have also analyzed the gaps and 
surface correlations of an S" = 1/2 chain with S = 1 im- 
purity end spins, which was predicted to be equivalent in 
the VBS picture. Indeed, we found that the exponent of 
the logarithmic decay of the correlations and the prcfac- 
tor of the logarithmically vanishing scaled first gap are 
universal. 

Higher gaps of the spectrum, in particular the second 
gap, scale as 1/L. The precise analysis of the second gap 
is numerically more demanding. Whereas we were able to 
check and more or less confirm some existing conjectures 
for the impurity models, we could not convincingly reach 
the thermodynamic limit for the S = 3/2 chain. Up to 
the lengths we could achieve in our DMRG calculation 
the anomalous dimension X2 shows a non-trivial value 
strictly below the predicted exponent a;2 = 1. More work 
is needed in the future on this issue. 

We close our paper with some remarks. First, we note 
that for not too long chains the numerical data can be 
easily misinterpreted to have an effective gap exponent, 
^ > 1, and a finite limiting end-to-end correlation func- 
tion, which is approached as a power-law. [For the anal- 
ogous inhomogeneous Ising quantum chain this happens 
for strong defect couplings, A > 1, see in Eqs. HA5|I and 
(|A4|) .] This type of behavior was suggested in an ear- 
lier numerical work;^ but data on longer chains gives 
evidence for the logarithmic dependence. 

Half- integer spin chains with S* > 3/2 are expected 
to have the same type of logarithmic end-to-end correla- 
tions as the S = 3/2 chain. In these cases the exponent 
of the logarithmic decay, 2d, is probably S dependent. 
We wonder if a field-theoretical approach, such as loga- 
rithmic conformal field theory can predict the values of 
these exponents. 

Our final remark concerns the effect of quenched dis- 
order. The quasi surface order in the S = 3/2 chain is 
found to be robust against weak quenched disorder. As 
studied in Rcf. [3^ the end-to-end correlations show the 
same qualitative behavior for not too strong randomness, 
and at the same time the critical properties of the sys- 
tem remains the same as in the non-random chain. For 
stronger, but still finite disorder the (quasi) surface or- 
der vanishes, and the critical behavior of the system is 
controlled by new types of disorder fixed pointsi^^*^ 
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APPENDIX A: SURFACE CRITICAL BEHAVIOR 
OF THE ISING QUANTUM CHAIN WITH AN 
EXTENDED DEFECT 

Here we consider an inhomogeneous Ising quantum 
chain, for which many properties such as the end-to-end 
correlation function and the excitation spectrum can be 
calculated exactlyj^^^i and these show similar character- 
istics as those of the 5 = 3/2 chain. These exact results 
are expected to serve as an analogy in the analysis. 

The inhomogeneous quantum Ising chain is defined by 
the Hamiltonian 



1 

Hr = -- 



x-i 



J2 J^'^i'^i'+i + J2 



with a smoothly varying coupling 



Ji = l 



A 


'L . 






— sm 




2 


TT 





(Al) 



(A2) 



in terms of the Pauli-spin operators, crf'^, at site 
This quantum Hamiltonian is related through a confor- 
mal transformation to a classical two-dimensional semi- 
infinite Ising model, the couplings of which at a distance 
I from the free surface deviate from the bulk value by an 
amount of A/ 21. 

Here wc list exact results about the end-to-end corre- 
lation function, Cl(1) = {afaf), as well as for the lowest 
gap ei, which can be calculated in the free-fermion rep- 
resentation, when the Hamiltonian is expressed as2i 



(A3) 



in terms of fermion creation, 77+, and annihilation, 77^, 
operators. The energy of the fermionic excitation, e^, is 
obtained from the solution of an L x L eigenvalue equa- 
tion: (A -I- B)(A — B)7/'q = EqV'qj whcrc the eigenvector, 
Tpq, contains information about the localized or delocal- 
ized nature of the excitation. 

The quantum Ising model in Eq. (|A1() has the same 
type of bulk quasi-long-range order as the homogeneous 
chain with ^ = 0, but the surface properties of the sys- 
tem are modified by the smoothly varying inhomogene- 
ity. These properties are qualitatively different for ^ < 1, 
A = \ and for A> I. 

The behavior of the end-to-end correlation function 
can be summarized as follows: 



Cl(1) 



const ^ 
const ln~^(L) 



if ^ < 1 
iiA=l 



(A4) 



1? + const L^--^ if y4 > 1, 
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i.e., it tends to zero as a power-law for ^ < 1, logarith- 
mically for A = I, and to a finite value with power-law 
corrections ior A > 1. 
The first gap scales as 



ei = < 



TT 1 



A 



2 

1 TT 1 

2lhrL 
, const L"-^ if ^ > 1. 



if A < 1 
if ^ = 1 



(A5) 



The gap vs exponent relation established by Eqs. (|12|l . 
(fTIjl and (|2SIl is satisfied for A <1 and A = l. For A> 1 
the first gap is determined by a mode localized to the 
surface. This latter becomes clear when the eigenvector 
of the lowest excitation, {~ytpi{l) is considered, which is 



given in the continuum approximation in terms of z 
ttI/L as 



1 sin^-^/^(z/2) 

Ll/2 C0S'4/2(z/2) 

tan^/2(^/2) 



if A < 1 

1 (A6) 
if A > 1. 



This eigenvector is localized at ^ = L for A > 1, loga- 
rithmically delocalized (quasi-localized) for A = 1 and 
delocalized for A < 1. Note that another surface excita- 
tion can be constructed, which is localized at / = 1, as 
dictated by symmetry. 
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